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Abstract 

Non-compact three-dimensional QED is studied by computer simulations to understand its 
chiral symmetry breaking features for Nf > 2, on lattice volumes up to 50^ and bare masses 
as low as rriQa = 0.0000625. We compute the chiral condensate, scalar and pseudoscalar 
susceptibilities, and the masses of scalar and pseudoscalar mesons. Finite volume effects 
and discretisation artifacts are carefully monitored. Our results reveal no decisive signal for 
chiral symmetry breaking for any Nf > 2. For Nf = 2 the dimensionless condensate can be 
bounded by < 5 x 10"^. We also present an exploratory study of the fractionally- 

charged Polyakov line. 
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1 Introduction 



Over the past few years QED3 has attracted a lot of attention, because of potential appli- 
cations to models of high superconductivity It is also an interesting and challenging 
model field theory in its own right, which is seen as an ideal laboratory to study effects 
such as dynamical mass generation which may be relevant to phenomena described by more 
complicated gauge field theories such as technicolor 0. In three dimensions the model is 
super-renormalisable, because the coupling g has mass dimension | and so provides the in- 
trinsic scale of the theory independently of any bare fermion mass mo- This also implies 
asymptotic freedom, since for processes with momentum transfer k ^ g"^ the theory is effec- 
tively non-interacting. 

Non-trivial behaviour may arise in the infra-red, however, as suggested by an expansion 
in 1/Nf, where Nj is the number of (four-component) fermion species In the large- A^j 
limit the photon propagator is modified by the vacuum polarisation diagram from l/k"^ to 
l/lk'^ + Y^fk]] in coordinate space the logarithmic Coulomb potential is thus modifed to 1/r 
for distances r ^ (g'^Nf)^^. Naively we deduce that the confining property of the Coulomb 
potential is screened by virtual // pairs. However, consider the energy required to excite a 
real // pair into the vacuum state [Q. The kinetic energy of the pair is positive and scales 
as by the uncertainty principle. For large r therefore, both kinetic and potential terms 
scale as and it becomes a delicate question which dominates. Next we note that the 
dimensionless interaction strength scales in this regime as Nj^, implying the possible exis- 
tence of a critical Nfc below which the negative potential term dominates as r 00. Since 
as r — »• the positive kinetic term must dominate the logarithmic Coulomb term, we deduce 
in this case the existence of an energy minimum at some non-zero r, implying the existence 
of stable // bound states in the ground state. This semi-classical argument hence suggests 
a non-vanishing condensate 7^ for Nf < Nfc- Since the condensate spontaneously 

breaks chiral symmetry, the fermions are imbued with a dynamically-generated mass. For 
asymptotically large r the massive fermions decouple, ceasing to screen the charge, and the 
logarithmically-confining Coulomb potential is restored. The picture therefore predicts two 
possible phases for QED3: for Nf < Nfc chiral symmetry is broken, so that the only massless 
states are Goldstone bosons, and electric charge is logarithmically confined; for Nf > Nfc 
the theory is conformal, consisting of massless fermions interacting via a long-range 1/r 
interaction. 

In this paper we address the question of whether chiral symmetry is spontaneously broken 
in QED3 for any value of Nf > 2, i.e. whether 7^ in the chiral limit bare fermion 

mass rriQ — > and the continuum (weak coupling) limit g 0. This question has received 
much study over the past decade, the principal approach being self-consistent solutions to 
truncated Schwinger- Dyson (SD) equation for the fermion propagator. Initial studies using 
the photon propagator derived from the leading order large- iVj expansion suggested that 
the answer is positive with Nfc — 3.2 Within this framework it was shown that the 
dynamically generated fermion mass is much less than the coupling contant, providing a 
natural hierarchy of mass scales. It is also possible to consider analytically continuing to a 
continuous number of flavors; the model in the limit Nf — > Nfc is then supposed to undergo 
an infinite-order phase transition . Other studies taking non-trivial vertex corrections into 
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account, however, predicted chiral symmetry breaking for arbitrary Nf [0. More recent 
studies which treat the vertex consistently in both numerator and denominator of the SD 
equations have found Nf^ < oo, with a value either in agreement with the original study 
P], or a slightly higher Nfc — 4.3 0. Other analytical approaches have also been applied. 
A recent renormalization group analysis predicts 3 < Nfc < 4 [jTO|. Finally there is an 
argument based on the inequality 

fiR < fuv. (1) 

where / is the absolute value of the thermodynamic free energy, and can be estimated in the 
relevant limits simply by counting massless degrees of freedom, either fermions or Goldstones 
as reviewed below. For QED3 this yields the prediction Nfc < | ||TT|. 



There have also been attempts to resolve the issue via numerical simulation of non- 
compact lattice QED3. The obvious advantage of a computer simulation is that one can study 
any Nf without any technical assumption concerning the convergence of expansion methods. 
Indeed, it is also possible to study the quenched {Nf = 0) limit, in which case a numerical 



study has shown that chiral symmetry appears to be broken |I2[, with a dimensionless 
condensate /?^('I'\E') ~ 5 x 10~^, where (3 = 1/g'^a is the dimensionless lattice coupling 
constant (a being the lattice spacing). For the unquenched theory, opinions have divided 
on whether Nfc is finite and ^ 3 @, |13|, or whether chiral symmetry is broken for all Nf 



14| . The most recent study with Nf = 2 on lattices up to 16 found chiral symmetry 



breaking with f3'^{'^'$) ~ 3 x 10 ^ |]T5|. This value is extracted, however, by simultaneous 



extrapolation to continuum, chiral, and thermodynamic limits, which is difficult to control 
since even with broken chiral symmetry there remain large finite volume effects due to the 
presence of a massless photon in the spectrum. For this reason it is unlikely that a numerical 
simulation will ever completely exclude chiral symmetry breaking for Nf > 0, particularly if 
the condensate is suppressed by a typical non-perturbative factor exp{—CNf) 0. 

In this paper we will present results of simulations on lattices considerably larger (up to 
50^) and bare fermion masses considerably smaller (moa > 0.0000625) than any previous 
study, which enables us to approach all three limits with greater control. Our goal is to 
place a bound on any possible chiral condensate using as few assumptions as possible on the 
nature of the required extrapolations. Our conclusion is that for Nf = 2, /3^(^\E') < 5 x 10~^, 
two orders of magnitude below the quenched value. For larger Nf the condensate is smaller 
still. Since we see no direct evidence for chiral symmetry breaking, our results are consistent 
with Nfc < 2, although as explained above an exponentially suppressed condensate is almost 
impossible to exclude. In the next section we briefly review the model, its symmetries, and 
the numerical methods used to simulate it. Our numerical results for the condensate appear 
in Section 3, in which in addition we present results for meson masses and susceptibilities, 
and Polyakov lines with both integer and fractional charges. All of these subsidiary quantities 
reveal evidence of strong finite volume effects which makes any conclusions drawn qualitative 
at best. Section 4 summarises our findings. 
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2 The Model 



The continuum Lagrangian density describing QED3 is given in Euclidean metric by 

£ = ^F^uF^^ + ^^D^-ff,^i + mo^i^i (2) 

where = 9^ — igAf^, and Ffj_^ is the field strength for the abelian gauge field A^. The 
fermions \E'j, z = 1, . . . , Nf, are four- component spinors acted on by hermitian Dirac matrices 
7^, implying that the mass term proportional to mo is invariant under the 2+1 dimensional 
parity transformation. 

Since the 7-matrices belong to a reducible representation of the Dirac algebra in three 
dimensions, the global symmetry of (|]) is larger than naively expected. The matrices 
{1, 73, 75, T3} where T3 = 27375 generate a global U(2) symmetry. To see this define \^ = ^t^. 
For mo = the U(2) symmetry is then 

^ ^ f/^ ^ ^ U e U(2). (3) 

The symmetry is broken, either explicitly by a bare mass or spontaneously by chiral con- 
densation (^'\E') 7^ 0, to a residual U(1)(8>U(1)t-3 where the first U(l) factor corresponds to 
global fermion number conservation, and the second to a symmetry 

^ e^»^3-q, _^ ^fe"'"'^^ (4) 

Chiral symmetry breaking for the theory with Nf fiavors would then result in the pattern 
Vi2Nf)^\]iNf)0ViNfU. 

The lattice action using staggered lattice fermion fields x, X is given by 

3 ^ 
S = ^ E ^M-(^)^M-(^) + EEXi(a;)M(x,x')x.(a;') (5) 

^ x,n<v 1=1 x,x' 

F^l/(^X^ = Ox^ ~\~ Ox+fi,u ^X + U^I^L ^xv 

The vectors x, x' consist of three integers [xq^xx^x-i) labelling the sites of an lattice. 
Since the gauge action is unbounded from above, @ defines the non-compact formulation 
of lattice QED3. The r7^(x) are Kawamoto-Smit phases {rii{x) = 1; r]2{x) = (— 1)"^^; tiq^x) = 
(^_iyi+x2'j designed to ensure relativistic covariance in the continuum limit. For the fermion 
fields antiperiodic boundary conditions are used in the timelike direction and periodic bound- 
ary conditions in the spatial directions. The phase factors in the fermion bilinear are defined 
by Uxfi = exp{i9xfj.), where 6*^.^ is the gauge potential. In terms of continuum quantities, 
Oxfi = cigAfj_{x) and the coupling /? = where a is the physical lattice spacing. We prefer 
the non-compact version of lattice QED3, because it resembles more closely the continuum 
model we wish to compare with. In particular, magnetic monopole excitations, which are 
known to have a decisive effect on the infra-red behaviour of the compact model are 
here strongly suppressed as /3 is made large. 
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The particle content of eq.(^ has been analysed in |T^, showing that in the long- 
wavelength limit flavors of staggered lattice fermion correspond to Nf = 2N of the 
four- component continuum flavors described by the action (Q). For lattice spacing a > 0, 
however, the global symmetries for mo = are only partially realised. In this case the 
symmetry is 

V — ¥ V P^°' V — > <=~*'^v 

AO Aoc- Ae 

Xe Xee'^ Xo e-^^Xo, (6) 

where Xo/e denotes the fleld on odd (i.e. e{x) = (— i^^i+^a+^s _ g^,^^ even sublattices 
respectively. Chiral symmetry breaking characterised by {^'^) = {xx) 7^ now has the 
pattern U(A^) U(A^) —>■ U(A^). Only in the continuum limit — >■ cxd is the full symmetry 
of the continuum theory expected to be restored. 

In the context of the current discussion the issue of whether the full U(2A^/)=U(4A^) 
flavor symmetry is restored in the continuum limit (3 —>■ oo is very important. Take for 
instance the considerations of ref. |Tl| which constrain Nf^ via the inequality (|1|): 

l + !(4iV^c)<l + 2iV|,. (7) 

Here the left hand side of the inequality counts the number of massless thermodynamic 
degrees of freedom in the UV limit, which consist of a single bose (the photon) and 4Nf 
fermi degrees of freedom, the factor | coming from the Fermi-Dirac distribution in 2-1-1 
dimensions. The right hand side counts both the photon and the 2Nj Goldstone modes 
expected from the breaking U(2A^/) U(A^/) ® U(A^/). Eq. (0) yields the bound Nf^ < f . 
If, however, chiral symmetry broke according to the staggered pattern U(A^)®U(A^) U(^), 
then the inequality becomes 

l + |(8iV,) < l + iV2, (8) 



implying Nfc = 2Nc < 12 ||l5l. An observation of broken chiral symmetry in the /? — oo 
limit in lattice QED3 simulations with Nf > 2 therefore implies either that the arguments 
leading to the constraint (P are incorrect, or that flavor symmetry restoration does not occur 
in this model. 

The numerical results presented here were obtained by simulating the action (|^) using 
a standard Hybrid Monte Carlo (HMC) algorithm. The form of (^) permits an even-odd 
partitioning, which means that a single flavor of staggered fermion can be simulated. The 
minimum number of continuum flavors which can be simulated with a local action is thus 
Nf = 2. We used the following methods to optimise the performance of the HMC algorithm. 
Firstly the effective coupling jS' used during integration of the equations of motion along a 
microcanonical trajectory was tuned so as to maximise the acceptance rate of the Monte 
Carlo procedure for a flxed microcanonical time-step dr. As the lattice size was increased dr 
had to be taken smaller and the optimal (3' approached (3. For example for the Nf = 2 case on 
a 16^ lattice and coupling P = 0.60 the choices dr = 0.025 and /?' = 0.6007 gave acceptance 
rate greater than 90% for all bare masses of interest. To maintain this acceptance rate on 
a 32^ lattice we used dr = 0.015 and /?' = 0.6004. Secondly the Monte Carlo procedure 
was optimised by choosing the trajectory length at random from a Poisson distribution 
with mean f = 1. This method of optimisation, which guarantees ergodicity, also decreases 
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autocorrelation times dramatically. Statistical errors in our measurements were calculated 
by jackknife blocking, which accounts for autocorrelations in a raw dataset. We found the 
errors remained stable as the number of blocks was varied from 10 to 50. 

3 Simulation Results 

3.1 Chiral Symmetry Breaking 

First let us make a statement of the problem we are addressing. Since a continuous symmetry 
is never broken spontaneously on a finite volume, we are required to work with mo 7^ and 
study the limiting behaviour as the explicit symmetry breaking is removed. Since for an 
asymptotically-free theory the UV behaviour is governed by the gaussian fixed point at the 
origin, then the continuum limit of the model lies in the limit (3 00, and all physical 
quantities should be expressible in terms of the scale set by the dimensionful coupling g. 
To compare simulation results taken at different couplings (lattice spacings), therefore, it is 
natural to work in terms of dimensionless variables such as jSmQ, L//3, or As the 

continuum limit is approached, data taken at different (3 should collapse onto a single curve 
when plotted in dimensionless units. Formally, chiral symmetry in QED3 is broken if 

lim lim lim B'^l^'^) ^ 

/3^oo/3mo^OL//3^oo 

In practice taking the limits in the required order is difficult, for several reasons: (i) 
The lattice itself distorts continuum physics considerably unless the lattice spacing a can be 
chosen small compared to the relevant physical wavelengths in the system - as we already 
mentioned in the previous section a non-zero lattice spacing breaks explicitly the continuum 
\J{2Nf) symmetry of the action; (ii) the size of the lattice must be large not just in 
dimensionless units but also compared to any dynamically generated correlations in the 
system; and (in) the chiral extrapolation tuq requires some theoretical prejudice. If the 
volume is large enough, mo lies in the linear regime of (^\l/(mo)) = (^^\l/(0)) + mo{^'^{0))' + 
. . . and linear extrapolations of the data work well. However, the chiral extrapolation is based 
on assumptions which are not necessarily valid when mo = 0. New physics at mo = may 
alter the situation. In our study despite the fact that in all our simulations we reached the 
"linear regime" we take a conservative stance by defining an upper bound for the condensate 
which is half its value at the smallest mo. 

In Fig. |I]we plot the dimensionless chiral condensate vs. the dimensionless bare 

mass /3mo for Nf = 2,4,8,16. The coupling (3 = 0.6 and the lattice volume is 16^. We 
generated approximately 1000 trajectories for each data point - statistical error bars are 
generally smaller than the size of the symbols. As Nf increases the chiral condensate de- 
creases (for mo > 0) because of increased screening of the interaction by dynamical fermions. 
As mo — >■ all the curves tend to pass smoothly through the origin. There is no sign of 
any discontinuous behaviour as Nf is altered. This motivates us to study in more detail the 
pattern of chiral symmetry breaking at small Nf on larger volumes near the chiral limit. 

To check whether lattice data are characteristic of the continuum limit we next plot in 
Fig. ^ the dimensionless chiral condensate /5^(^^\I') vs. the dimensionless fermion bare mass 
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(9) 



(3mQ for A^j = 2 at different values of the coupling (3 = 0.45, 0.60, 0.75, 0.90. In order to 
disentangle the lattice discretisation effects from the finite size effects we keep the volume 
in physical units (L//3)^ constant. It can be easily inferred from the graph that at strong 
coupling (/3 = 0.45) discretisation effects are significant whereas for (3 > 0.60 the lattice 
artifacts become small since the data almost fall on the same line within the resolution of 
our analysis. 

Having clarified the issue of the continuum limit we next investigate the thermodynamic 
limit, simulating with Nf = 2 aX (3 = 0.6. In Fig. ^ we present our results for (3'^{^'^) 
vs. /3mo on different lattice sizes: 8^, 16^, 24^, 32^ and 48^. Note that the volumes and 
masses explored here are considerably closer to the thermodynamic and chiral limits than 



those of the study of |T5|, for which L/(3 < 32 and (3mo > 0.005. We generated 300 - 800 
configurations for each data point. This enables us to expose finite size effects at different 
values of the fermion bare mass. Finite size effects become small (but still significant) for 
L > 24 {L/f3 > 40), which is consistent with the quenched results of |]T2|; however as L 
increases all the lines continue to pass smoothly through the origin, suggesting that chiral 
symmetry remains unbroken. 

As well as the chiral condensate, we have also studied mesonic correlation functions. In 
many cases these quantities yield useful information about the nature of the ground state 
model, although they are also more prone to finite volume effects than the simple order 
parameter. Fig. ^ shows results from Nf = 2, (3 = 0.6 on lattice size 16^ and 32^ for the 
ratio of longitudunal to transverse susceptibilities R = Xi/Xt, where xi,t are the integrated 
propagators in scalar and pseudoscalar meson channels respectively: 

xi = J2(xx{o)xx{x)); xt = J2(x£x{o)x£x{x)). (10) 

X X 

The longitudinal susceptibility xi has contributions from diagrams with both connected and 
disconnected fermion lines; we checked that the connected contribution is dominant through- 
out our parameter space. These quantities are much noisier than the chiral condensate. The 
transverse susceptibility is most conveniently estimated via the Ward identity Xt = {^'^) /mQ. 
In the chiral limit we expect 

^. ^ I 0, chiral symmetry broken; 
mo^o 1 1, chiral symmetry unbroken. 

As we can see from the figure, at large values of mo the finite size effects are small, but 
become significantly larger at intermediate values of mo. This implies that finite volume 
effects are not fully under control; the fact that R from both 16^ and 32^ lattices appears to 
converge to 1 as mo — >■ decreases does not necessarily imply chiral symmetry is unbroken. 
It could be attributed to finite volume effects. Nonetheless, there is no indication of any 
range of mo over which dR/drriQ is positive, required if chiral sjTiimetry is broken. In Fig. ^ 
we plot R vs. mo for Nf = 2, 4. This result is compatible with the data presented in Fig. |I], 
i.e. that the explicit chiral symmetry breaking at nonzero mo, indicated by the departure of 
R from unity, is stronger as Nf decreases. 

In Fig. § we present the results for the pion mass M^^ vs. fermion bare mass mo for 
Nf = 2, (3 = 0.6 and lattice sizes 16'^, 32^, 48^. If the model is chirally symmetric we expect 
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cx mo; otherwise from conventional PCAC arguments we expect M^. oc ^Jrn^. As in the 
case of i?, suffers from very strong finite volume effects and therefore decisive conclusions 
cannot be drawn. In Fig.|^ we plot the pion and sigma masses vs. mg for N ^ = 2, /? = 0.6 and 
lattice size 32'^. The masses tend to become equal as mo — > 0, and show no sign of tending 
to zero in the chiral limit. It should be noted in this regard that there are SD calculations 



predicting the absence of light scalar excitations for Nj > Nf^ ||T8|. However, as already 
stressed no definitive conclusions can be drawn until finite size effects are under control. 

Since our best chance of controlling the extrapolation to infinite volume appears to come 
from measurements of the order parameter, where signals are less noisy and finite volume 
effects relatively small, we now return to this and discuss results from A'^^ = 2 simulations 
at /5 = 0.75, closer to the continuum limit than our previous simulations. These datasets 
result from between 300 and 500 HMC trajectories. We studied the behaviour of the chiral 
condensate in detail in the large volume and small mg limits. Fig. ^ shows results for 
^2|/i|rv|/^ vs. l3mo extracted from simulations on system sizes 10'^, 20^, 30^, 40^ and 50^. These 
simulations were performed very close to the chiral limit (mo < 0.005). Although we don't 
extrapolate our data to the chiral limit it is clear that all the data tend to pass smoothly 
through the origin. A comparison of the 40^ and 50^ points shows that finite volume effects 
are under control, at least for Puiq > 0.0007. Fig. |^ shows the same data zoomed to the 
vicinity of the origin. We can see that even when mo is very small, the data tend to pass 
smoothly through the origin: there is no indication of any tendency for the extrapolated 
intercept to be non-zero signalling chiral symmetry breaking. With the conservative criterion 
adopted above we conclude that for Nf = 2, /3^(\E'\E') < 5 x 10~^, which is a strong indication 
that chiral symmetry in QED3 may be restored for Nf > 2. 

In Figs. H and |] we plotted the condensate vs. bare mass keeping the lattice volume 
constant along each line. Hitherto we have assumed that the physical volume of the system 
is proportional to {L/j3)^. If dynamical mass generation occurs by some means, however, so 
that a correlation length ,^ < 00 develops, then the correct measure of the physical volume 
is the ratio {L/^Y (this is merely the statement that finite volume effects are sensitive to 
which phase the theory resides in). This quantity varies along the constant lattice volume 
lines of Figs. because ^ is in principle a function of both mo and (3. We have attempted 
to study finite volume effects taking this possibility into account: since however we don't 
have an accurate estimate of a correlation length ^ either from the sigma mass, because of 
strong finite volume effects, or from the fermion mass, because the corresponding propagator 
is not gauge invariant, we are forced to estimate the behaviour of ^(/5, mo) directly from the 
condensate data. Since the precise relation between ^ and {^'^) is unknown, we proceed by 
observing that in all of our datasets for fixed (3 and mo > 0.001 it is approximately true that 
oc mo and assuming that ^ oc (^'\I')^" oc mg " for some power a = {d—l—ri^q,)~^ where 
rj^q, is the anomalous scaling dimension of the composite operator We then cover the 
range of possibilities by evaluating the condensate keeping LttLq constant, with the choices 
a = \ expected for a QCD-like theory in which r^^,j, is perturbatively small; and a = 1, 
expected for a theory such as the NJL model in which r]^^ ^ d — 2 and the condensate is 
directly proportional to the dynamically-generated scale. A value predicted for QED3 with 
Nf < Nfc by the SD approach is a = | corresponding to 77^,^ = | 0- 

The results for /? = 0.75 and mo > are shown in Figs. ^ and |ll] respectively. In the 
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case of the constant Ly/rn^ normalisation the lattice sizes for the three curves shown are 
[503,363,303,263,223], [403,283,243,203,183] and [303,223,183,163,143]. For constant Lmo 
the sequences are [503, 263, 183, 143, ^gS]^ [4o3, 203, 143, iq3^ gS] [ggs^ ^gs^ ^03, 83, 63]. By 
the very nature of the way the curves are produced an extrapolation to a non-trivial chiral 
limit is impossible, since eventually the assumption that oc mo must break down; we 

note, however, that both sets of curves have negative curvature. It is also apparent that 
the curves with constant L^itiq appear to approach the infinite volume limit uniformly over 
a wide range of mo, whereas those keeping Lmo constant diverge as mo increases. This is 
tentative evidence against a large value for 77^,^. 



3.2 Deconfinement of Fractional Charge 

Finally we consider a different issue, namely the confining properties of the theory. Since the 
gauge degrees of freedom are non-compact, matter fields and/or test charges may in principle 
be defined with any value of electric charge g which need not be commensurate with the 
dynamical fermion charge g. If the latter is the case, then such fractional charges should be 
confined by the Coulomb potential independent of the phase of the model, since long-range 
forces cannot in this case be screened by virtual // pairs. It has been pointed out, however, 
that when the chiral symmetry of QED3 is broken, either explicitly or spontaneously, then 
there may be a transition to a phase where fractional charge is deconfined at non-zero 
temperature [ITU . 



In brief, on a system with finite extent in the Euclidean time direction, the action 
(@) has a global Z symmetry under non-periodic gauge transformations; in a static gauge 
this corresponds to invariance under 

nT 

Ao{x) ^ Ao{x) + 27r— , n e Z. (12) 

It is readily seen that worldlines of all integer-charged matter fields are invariant under 
(p!^). For fractionally- charged particles, however, this is no longer the case. The analysis 



of [|T9l derives an effective potential V{Ao) which is periodic under (|l^) ; as T increases the 
barrier between non-equivalent minima of V grows until the system becomes trapped near 
one such minimum, breaking the Z symmetry, with the fractionally- charged Polyakov line 
Ilg{x) = expig Jq AQ{x)dxo acting as a gauge invariant order parameter. By analogy with 
Yang-Mills theory the symmetry broken phase is identified as one in which fractional charge 
is deconfined. 

In the simulations of this paper we have not examined systems with the Euclidean time 
direction distinguished; however, we have found it interesting to check whether the Z sym- 
metry is broken as a means of monitoring finite volume effects. We have chosen to examine 
charges which are rational fractions of the fundamental charge: g = g/n, and define the 
Polyakov loop Iln{x) by 

nr,ix) = l[exp(-eoix,t)). (13) 



t=l 



n 



The original Z symmetry (|T2D translates into a Z„, symmetry on n„; spontaneous breaking 
is therefore signalled by a non-zero expectation of Re(L~3 ^^n„(x))". In Fig. ^ we plot 
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distributions of this quantity for n = 1,2,4 and 8 on a sequence of lattice volumes. On all 
systems the n = 1 data is sharply peaked about zero, indicating that there is no deconfine- 
ment of the fundamental integer charge|^. On 8^, however, the distributions for all n > 2 
are skewed in the positive direction, indicating deconfinement. As the volume is increased, 
permitting fluctuations of the flelds to develop over longer distances, the deconflnement sig- 
nal is washed out until on 32'^ only charges with g = g/8 remain deconflned. A possible 
explanation of why smaller charges are harder to conflne in a flnite volume is that g"^ itself 
provides the scale in the logarithmic potential. 



4 Summary and Outlook 

In our study of QED3 with Nf > 2 we attempted to establish whether chiral symmetry 
is broken or not by studying various observables close to the continuum limit (7 — > 0, on 
different volumes in order to detect and control flnite size effects and near the chiral limit 
mo — >■ 0. As expected, susceptibilities and meson masses suffer from strong flnite size effects 
and by themselves they do not allow us to reach a deflnitive conculsion. However, for the 
chiral condensate the continuum, thermodynamic and chiral limits are under better control. 
On the volumes we have been able to study we have seen no evidence for chiral symmetry 
breaking for any Nf > 2. 

Our upper bound for the condensate in the Nf = 2 case is < 5 x 10~^, to be 

compared with the estimate for the quenched model /3^(\E'\E') ^ 5 x 10"'^ @, |l2l. In addition. 



for all Nf > 2 the curves of /?^(\E'\E') vs. /3mo tend to pass smoothly through the origin 
with no sign of any discontinuous behaviour in Nf. This is evidence against the suggestion 
that the IR flxed point behaviour of QED3 for Nf < Nf^ coincides with the UV behaviour 
of the three dimensional Thirring model as suggested by both the large- A^j expansion and 
SD studies . For the Thirring model the critical Nf^ below which a non-trivial UV flxed 
point exists has been estimated by Monte Carlo simulation to be Nfc — 5 |20, 21|. The 



two models appear to lie in different universality classes, possibly because for Nf < Nfc the 
Thirring model has no massless degree of freedom corresponding to the photon. 

We are currently simulating the model with Nf = 1 to check whether in this case chiral 
symmetry is broken, as predicted by |Tl|. It will also be interesting in this respect to study 



whether the continuum \J{2Nf) symmetry is restored; this may entail the use of improved 
fermion actions |^ in future simulations. We also plan to extend our work to study in detail 



the deconflnement phase transition of fractional electric charges at T > 0. 



^More precisely, it suggests the free energy of an isolated fundamental charge diverges in the thermody- 
namic limit, which could also occur in the conformal phase. 



10 



Acknowledgements 

SJH and CGS were supported by a Leverhulme Trust grant, and JBK in part by NSF grant 
PHY-0102409. The computer simulations were done on the Cray SVl's at NERSC, the 
IBM-SP at NPACI, and on the SGI Origin 2000 at the University of Wales Swansea. We 
have enjoyed discussions with Ian Aitchison, Nick Mavromatos, Sarben Sarkar and Rohana 
Wijewardhana. 

References 

[1] N. Dorey and N.E. Mavromatos, Phys. Lett. B250 (1990) 907; Nucl. Phys. B386 (1992) 
614; 

K. Farakos and N.E. Mavromatos, Mod. Phys. Lett. A13 (1998) 1019, 

I.J.R Aitchison, N. Dorey, M. Klein-Kreisler and N.E. Mavromatos, Phys. Lett. B294 

(1992) 91; 

I.J.R Aitchison, Z. Phys. C67 (1995) 303. 

[2] B. Holdom, Phys. Lett. B150 (1985) 301; 

K. Yamawaki, M. Bando and K. Matumoto, Phys. Rev. Lett. 56 (1986) 1335; 

T. Appelquist, D. Karabah and L.C.R. Wijewardhana, Phys. Rev. Lett. 57 (1986) 957. 

[3] R.D. Pisarski, Phys. Rev. D29 (1984) 2423. 

[4] E. Dagotto, A. Kocic and J.B. Kogut, Nucl. Phys. B334 (1990) 279. 

[5] T.W. Appelquist, M. Bowick, D. Karabali and L.C.R. Wijewardhana, Phys. Rev. D33 
(1986) 3704; 

T.W. Appelquist, D. Nash and L.C.R. Wijewardhana, Phys. Rev. Lett. 60 (1988) 2575. 

[6] V.A. Miransky and K. Yamawaki, Phys. Rev. D55 (1997) 5051; 

V.P. Gusynin, V.A. Miranskii and A.V. Shpagin, Phys. Rev. D58 (1998) 085023. 

[7] M.R. Pennington and S.R. Webb, Brookhaven preprint BNL-40886; 
M.R. Pennington and D. Walsh, Phys. Lett. B253 (1991) 246. 

[8] P. Maris, Phys. Rev. D54 (1996) 4049. 

[9] D. Nash, Phys. Rev. Lett. 62 (1989) 3024; 

K.-I. Kondo, T. Ebihara, T. lizuka and E. Tanaka, Nucl. Phys. B434 (1995) 85; 
I.J.R. Aitchison, N.E. Mavromatos and D.O. McNeiU, Phys. Lett. B402 (1997) 154. 

[10] K.-I. Kubota and H. Terao, Prog. Theor. Phys. 105 (2001) 809. 

[11] T.W. Appelquist, A.G. Cohen and M. Schmaltz, Phys. Rev. D60 (1999) 045003. 

[12] S.J. Hands and J.B. Kogut, Nucl. Phys. B335 (1990) 455. 



11 



[13] E. Dagotto, A. Kocic and J.B. Kogut, Phys. Rev. Lett. 62 (1989) 1083; 
J.B. Kogut and J.-F. Lagae, Nucl. Phys. B(Proc. Suppl.)30 (1993) 737. 

[14] V. Azcoiti and X.-Q. Luo, Mod. Phys. Lett. A8 (1993) 3635. 

[15] J. Alexandre, K. Farakos, S.J. Hands, G. Koutsoumbas and S.E. Morrison, Phys. Rev. 
D64 (2001) 034502. 

[16] A.M. Polyakov, Nucl. Phys. B120 (1977) 429; 

T. Banks, R. Myerson and J.B. Kogut, Nucl. Phys. B129 (1977) 493; 
R.J. Wensley and J.D. Stack, Phys. Rev. Lett. 63 (1989) 1764; 

M.N. Chernodub, E.M. Ilgenfritz and A. Schiller, Phys. Rev. Lett. 88 (2002) 231601. 
[17] C. Burden and A.N. Burkitt, Europhys. Lett. 3 (1987) 545. 

[18] T.W. Appelquist, J. Terning and L.C.R. Wijewardhana, Phys. Rev. Lett. 75 (1995) 
2081. 

[19] G. Grignani, G. Semcnoff and P. Sodano, Phys. Rev. D53 (1996) 7157; 
LJ.R. Aitchison and CD. Fosco, Nucl. Phys. B578 (2000) 199. 

[20] L. Del Debbio, S.J. Hands and J.C. Mehegan, Nucl. Phys. B502 (1997) 269. 

[21] L. Del Debbio and S.J. Hands, Nucl. Phys. B552 (1999) 339; 
S.J. Hands and B. Lucini, Phys. Lett. B461 (1999) 263. 

[22] C. Bernard et al (the MILC collaboration), Phys. Rev. D58 (1998) 014503. 



12 




0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 

l3mo 

Figure 1: Dimensionless condensate vs. dimensionless bare mass l3mQ for Nj = 

2, 4, 8, 16, (3 = 0.6 on a 16^ lattice. 




Figure 2: vs. (3mo for Nf = 2 at different values of the coupling (3 and constant 

physical volume (L//3)^. 
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Figure 3: vs. pmo for Nf = 2, p = 0.6 and lattice sizes 8^, 16^, 24^, 32^, 48^. 
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Figure 5: i? vs. mo for Nf = 2,A, (3 = 0.6 on 32^ lattices. 
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Figure 6: Pion mass vs. mo for iVy = 2, /3 = 0.6 on lattice sizes 16 , 32 and 48 . 
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Figure 8: Condensate vs. bare mass for Nj = 2, P = 0.75 and lattice sizes 10^,20^,30^,40^ 
and 50^. 
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Figure 9: Same as Fig.|3 but zoomed near the origin. 
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Figure 10: Condensate vs. fermion bare mass for Nf = 2 and /3 = 0.75 with L^frn^ kept 
approximately constant along each curve. The points at the smallest value of (5m^ were 
extracted from simulations on 50^, 40^ and 30^ lattices for each of the three different curves. 
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Figure 11: Condensate vs. fermion bare mass for N ^ = 2 and (3 = 0.75 with Lmo kept 
approximately constant along each curve. 
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Figure 12: Histograms of the distribution of Re(n„)"^ from simulations with Nf = 2, f3 = 0.6 
on 8^, 16^, 24^ and 32^ lattices, from respectively 10000, 2000, 1400 and 800 configurations. 
The squared shading denotes n = 1, cross-hatched to north-east n = 2, to north-west n = 4, 
and scalloped n = 8. 
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